SCIENTIFIC 

REPpRTS 



r>3 




OPEN 



SUBJECT AREAS: 

ELECTRONIC PROPERTIES 
AND MATERIALS 

TOPOLOGICAL INSULATORS 



Received 
27 May 2014 

Accepted 
30 July 2014 

Published 
20 August 2014 



Correspondence and 
requests for materials 
should be addressed to 
XJ.Z. (XJZhou@aphy. 

iphy.ac.cn) 



Evidence of Topological Surface State in 
Three-Dimensional Dirac Semimetal 
Cd 3 As 2 

Hemian Yi 1 , Zhijun Wang 1 , Chaoyu Chen 1 , Youguo Shi 1 , Ya Feng 1 , Aiji Liang 1 , Zhuojin Xie 1 , Shaolong He 1 , 
Junfeng He 1 , Yingying Peng 1 , Xu Liu 1 , Yan Liu 1 , Lin Zhao 1 , Guodong Liu 1 , Xiaoli Dong 1 , Jun Zhang 1 , 
M. Nakatake 3 , M. Arita 3 , K. Shimada 3 , H. Namatame 3 , M. Taniguchi 3 , Zuyan Xu 4 , Chuangtian Chen 4 , 
Xi Dai 1 , Zhong Fang 1 &X. J. Zhou 1 - 2 

1 Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 1 001 90, 
China, 2 Collaborative Innovation Center of Quantum Matter, Beijing, China, 3 Hiroshima Synchrotron Radiation Center, Hiroshima 
University, Hiroshima 739-8526, Japan, 4 Technical Institute of Physics and Chemistry, Chinese Academy of Sciences, Beijing 
1001 90, China. 

The three-dimensional topological semimetals represent a new quantum state of matter. Distinct from the 
surface state in the topological insulators that exhibits linear dispersion in two-dimensional momentum 
plane, the three-dimensional semimetals host bulk band dispersions linearly along all directions. In addition 
to the gapless points in the bulk, the three-dimensional Weyl/Dirac semimetals are also characterized by 
"topologically protected" surface state with Fermi arcs on their surface. While Cd 3 As 2 is proposed to be a 
viable candidate of a Dirac semimetal, more investigations are necessary to pin down its nature. In 
particular, the topological surface state, the hallmark of the three-dimensional semimetal, has not been 
observed in Cd 3 As 2 . Here we report the electronic structure of Cd 3 As 2 investigated by angle-resolved 
photoemission measurements on the (112) crystal surface and detailed band structure calculations. The 
measured Fermi surface and band structure show a good agreement with the band structure calculations 
with two bulk Dirac-like bands approaching the Fermi level and forming Dirac points near the Brillouin 
zone center. Moreover, the topological surface state with a linear dispersion approaching the Fermi level is 
identified for the first time. These results provide experimental indications on the nature of topologically 
non-trivial three-dimensional Dirac cones in Cd 3 As 2 . 



The Dirac materials, such as graphene and topological insulators, exhibit novel quantum phenomena and 
properties with their low energy electron behaviors described by the relativistic Dirac equation 1 ' 2 . The 
massless Dirac fermions on the surface of graphene 3 " 5 and the three-dimensional topological insulators 6 ' 7 
show cone-shaped band dispersion with the Dirac points lying at discrete points in two-dimensional (2D) 
momentum plane. Recently, much attention has been attracted to the Weyl/Dirac semimetals which represent 
a new topological state of quantum matter that is different from the topological insulators 8 " 17 . These so-called 
three-dimensional (3D) Weyl/Dirac semimetals have inverted band structures in which the bulk valence band 
and conduction band contact at discrete points (Dirac points) in three-dimensional momentum space. They can 
be viewed as the 3D analogue of graphene and the surface state of the topological insulators that host Dirac band 
dispersions linearly along all directions 1418 . In addition to the gapless points (Weyl/Dirac nodes) in the bulk that 
are "topologically protected" by the behavior of the band structure, the Weyl/Dirac semimetal is also character- 
ized by "topologically protected" surface state with Fermi arcs on the surface 910 . The Weyl/Dirac semimetal is 
expected to show various new quantum phenomena and unusual properties like magnetic monopole in the bulk, 
giant diamagnetism, quantum spin Hall effect or quantum anomalous Hall effect in its quantum well struc- 
ture 8 ' 91519 ' 20 . In the case with both the time reversal and inversion symmetries, one may expect a 3D Dirac 
semimetal state described as four- component Dirac fermions which can be viewed as two copies of distinct 
Weyl fermions 21 . To get Weyl semimetal, either the time reversal symmetry or the inversion symmetry has to be 
broken. 

While a number of materials and systems have been proposed for realizing the Weyl/Dirac semimetal, few has 
been materialized so far 910141518 ' 22 ' 23 . Among the proposals, Cd 3 As 2 has been paid a special attention for a couple 
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of reasons 18 . First, band structure calculations predict that Cd 3 As 2 is 
a symmetry-protected topological semimetal, with a single pair of 3D 
Dirac points in the bulk and nontrivial Fermi arcs on the surface 18 . 
The pseudo-relativistic physics of dispersions results in 3D Dirac 
points that are protected by crystal symmetry and robust against 
perturbations. No gap can be opened by the spin- orbit coupling, in 
contrast to Dirac materials like AMnBi 2 (A = Sr,Ca) 24,25 . Second, 
Cd 3 As 2 is chemically stable in air which makes it more convenient 
for property investigations and potential applications when com- 
pared to other candidates 15,26,27 . Third, similar to graphene, Cd 3 As 2 
has particularly high carrier mobility up to 15000 cm 2 V~ l s~ l at room 
temperature and 80000 cm 2 V~ l s~ l at 4 K 28,29 , which makes it ideal 
for revealing novel quantum phenomena from transport measure- 
ments as well as designing electronic devices. It provides a promising 
platform to realize quantum spin Hall effect and sizable linear 
quantum magnetoresistance even up to room temperature 18 . While 
the indications on the 3D Dirac fermion nature of Cd 3 As 2 are accu- 
mulating, with some angle- resolved photoemission results reported 
for the Cd 3 As 2 (001) surface 30,31 and (111) surface 32 , more experi- 
mental evidence and theoretical investigation are necessary to pin 
down the issue. In particular, the topological surface state, the hall- 
mark of the 3D semimetal, has not been identified in Cd 3 As 2 . If the 
3D Dirac fermion nature of Cd 3 As 2 can be fully proven, it will open 
an opportunity to expand further into a Weyl semimetal when the 
time-reversal or inversion symmetry is broken. 

Angle- resolved photoemission spectroscopy (ARPES) is a power- 
ful tool which can directly reveal the presence of the Dirac cone and 
its 3D character. Here we report our ARPES study and band structure 
calculations on the (112) crystal surface of Cd 3 As 2 . Linear-dispersive 
energy bands are observed in the low energy region that are assoc- 
iated with the bulk Dirac nodes. Such linear- dispersive behavior of 
the bulk bands is further examined along the out-of-plane direction 
by varying the photon energy of the light source. In particular, we 
have identified for the first time the topological surface state in the 
gapless bulk Dirac nodes of Cd 3 As 2 . These results are consistent with 
the band structure calculations and provide experimental indications 
of topologically non-trivial 3D Dirac cones in Cd 3 As 2 . 

Figure 1 shows the crystal structure, Brillouin zone and X-ray 
diffraction patterns of Ca 3 As 2 (see Supplementary for experimental 
details). The crystal structure of Cd 3 As 2 is body-centered tetragonal 
with an 14 led symmetry; it is related to the tetragonally-distorted 
anti-fluorite structure with 1/4 Cd site of ordered vacancy 33 . As seen 
from the small cubic sub-cell (bottom of Fig. la), the As ions are 
approximately close-packed and the Cd ions are tetrahedrally coor- 
dinated. The overall crystal structure of Cd 3 As 2 (top of Fig. la) can be 
viewed as a 2 X 2 X 4 superstructure stacked from the small cubic 
sub-cell. Fig. lb shows the reduced Brillouin zone and the projected 
(112) surface Brillouin zone for the small unit cell of Cd 3 As 2 that is 
used in our band structure calculations. Powder X-ray diffraction 
pattern (Fig. lc) on the sample from grinding the grown Cd 3 As 2 
single crystals provides an accurate determination of the lattice para- 
meters a = b = 12.6467 A and c = 25.4428 A. Fig. Id displays the 
XRD pattern of the Cd 3 As 2 single crystal. The observed peaks can be 
indexed into (n n 2n) (with n being integers) reflections which indi- 
cate that the naturally cleaved surface is ( 1 12). It actually corresponds 
to the close-packed As planes in the small cubic sub-cell with a 
distance of 3.68 A between two adjacent planes. 

Figure 2 shows constant energy contours of the spectral weight 
distribution for the cleaved (112) surface of Cd 3 As 2 crystal over a 
large momentum space (Fig. 2a), and their comparison with band 
structure calculations (Fig. 2b). Here the plotted hexagonal Brillouin 
zone is the projected (112) surface Brillouin zone (Fig. lb). The 
spectral weight repetition over three zones reaffirms that the mea- 
sured crystal surface is (112). Right at the Fermi level, the measured 
spectral weight is point-like that concentrates at the Brilllouin zone 
centers (r points) (Fig. 2al). With increasing binding energy, the 



spectral weight distribution gets expanded and at high binding ener- 
gies like 450 meV (Fig. 2a4) and 600 meV (Fig. 2a5), two fine fea- 
tures, the inner part and outer part, can be discerned as marked by 
blue dashed lines. The inner part appears as a small circle. The outer 
part changes from a point-like Fermi surface to a ring-like circle and 
eventually to a hexagon -like shape at high binding energy. From the 
band structure calculations, there are only Fermi points at the Fermi 
level (Fig. 2bl). With increasing binding energy, the point gradually 
becomes hexagon -shaped at 150 meV, and warped hexagon-shape at 
high binding energies. The spectral weight evolution with increasing 
binding energy is reminiscent with many materials with Dirac cones 
such as graphene and topological insulators. The hexagon -shape 
indicates that the related Dirac cone is anisotropic, a situation similar 
to that found in Bi 2 Te 3 34,35 . Although the calculated fine- details of the 
electronic structure are not fully resolved in our measurements, our 
results agree well with the band structure calculations. Overall, we 
find that the periodicity of the ARPES signal can be well described by 
using the small sub-cell (Fig. la); the effect of the (2X2X4) super- 
cell that is supposed to cause band folding is rather weak. 

Figure 3 shows band structures of Cd 3 As 2 (112) surface measured 
along two high- symmetry directions. There is a good agreement 
between the measured band structures (Figs. 3b and 3c) and the band 
structure calculations (Fig. 3d). Along the T-K cut, from the calcula- 
tion, there are four main bulk bands, labeled as BV1, BV2, BV3 and 
BV4, within the energy range shown (Fig. 3d). Along the T-M cut, 
there are three bulk bands, labeled as BV5, BV6 and BV7(Fig. 3d). As 
seen from the measured band structures (Figs. 3b and 3c) along T-K 
(cut B, green line in Fig. 3a) and T-M (cut A, red line in Fig. 3a), all 
these seven bulk bands from the calculations can find good corres- 
pondence. It is clear that, within the similar energy range, we have 
resolved more bands than reported before on the (112) surface of 
Cd 3 As 2 32 . The good agreement between the measured and calculated 
band width also indicates that the band renormalization effect in 
Cd 3 As 2 is weak. Furthermore, there is a surface band from the cal- 
culated results (labeled SV1 in Fig. 3d) that appears near the zone 
center (r point). The surface state only appears in surface band 
calculation when taking into account the spin-orbit coupling as well 
as the inverted bulk band structure near the Fermi level. As seen from 
Figs. 3b and 3c, there is a band in the measured data labeled as SV1 
that is in good agreement with the calculated surface band. We note 
that there is some band intensity asymmetry along the equivalent 
momentum cut locations. For example, as seen in Fig. 3b, along T x - 
M, three bands are clearly seen but along the equivalent r 21 -M, the 
BV5 band gets stronger and the BV6 band becomes invisible. For the 
SV1 band, it appears clearly only on one side of T-M (Fig. 3b) and T- 
K (Fig. 3c) measurements. These intensity variations can possibly be 
attributed to the photoemission matrix element effect. From the 
measurements, there are three bands, two bulk bands(BV2 and 
BV5) and one surface band (SV1), that are Dirac-cone-like 
approaching the Fermi level. The good agreement between the mea- 
surements and band structure calculations provide an unambiguous 
identification of each band in Cd 3 As 2 that lays a foundation for our 
further study below. 

The nature of the surface band SV1 can be further examined by 
varying the photon energy during the ARPES measurements to 
detect band structures at different K z planes. As seen in Fig. 4a, at 
different photon energies, the SV1 band (marked as dashed blue 
lines) shows little change within the experimental uncertainty. This 
again provides an additional strong evidence on the nature of the SV1 
surface band. We also find that this surface state is rather robust as it 
does not show obvious change during long time ARPES measure- 
ment even at a high temperature (—300 K). Interestingly, the SV1 
surface band lies within the gap region formed by the BV2 and BV3 
bands along the T-K and BV5 and BV6 bands along T-M. This is 
distinct from the topological surface state in 3D topological insula- 
tors such as Bi 2 Se 3 36 and Bi 2 Te 3 6 . This surface state possesses a Fermi 
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Figure 1 | Crystal structure, Brillouin zone and the cleaved surface of Cd 3 As 2 . (a). Crystal structure of Cd 3 As 2 with a body-centered tetragonal 
structure (Space group: 14 led). Such a unit cell can be viewed as a re-constructed 2X2X4 superstructure from a small sub-cell expanded in the bottom- 
left. The As ions form a face-centered cubic while the Cd ions fill 3/4 of the 8 tetragonal sites formed by the As ions. (b). The bulk Brillouin zone and 
projected (112) surface Brillouin zone used in the band structure calculations which corresponds to the smallest unit cell marked in the sub-cell as red 
arrows in (a), (c). The X-ray diffraction pattern of Cd 3 As 2 powder taken with Cu Ka radiation = 1.5418 A). The measured spectrum (red squares) is 
compared with the calculated result (green line) from the structure refinement. It confirms the crystal structure with a space group 14 led; the lattice 
constants obtained are a = b = 12.6467 A and c = 25.4428 A. (d). X-ray diffraction pattern of a naturally cleaved surface of the Cd 3 As 2 single crystal that 
indicates the (112) surface. 



velocity as high as 7.5 eV which has nearly the same value along G-K 
and G-M directions, suggesting an isotropic nature of the SV1 sur- 
face state. These results are consistent with the existence of a two- 
dimensional Dirac cone derived from the SV1 surface band. The 
revelation of the topological surface state is in agreement with the 
expectation from theoretical calculations on the non -trivial bulk 
band topology of Cd 3 As 2 18 . 

The existence of possible 3D Dirac fermion in Cd 3 As 2 , which is 
supposed to disperse linearly along all three directions in momentum 
space, can be checked by the two-dimensional energy dispersions in 
the K x -K 7 basal plane by varying K y values (Fig. 4b) and the band 
dispersion along K z by varying the photon energy (Fig. 4a). As seen 
from Fig. 4b, although the signal is weak, there is an overall trend 
that, with the momentum cut moving along K y within the K x -K y 
basal plane (cuts Bl to B5 in Fig. 4c), the cross point of the BV1 bulk 
band approaches the Fermi level at K y = 0 and sinks below the Fermi 
level when the momentum cut moves away from the T-K direction. 



The dispersion can be fitted with a linear line at the Brillouin zone 
center (cut 3). Away from the Dirac point, the dispersions can be 
fitted with hyperbola lineshape, similar to the Dirac cone formed on 
the surface of topological insulators 6,7 . The measured momentum 
dependence is consistent with the calculated result (Fig. 4d). These 
results are consistent with the presence of two-dimensional Dirac 
fermion from the BV1 band in the K^-K^ plane. On the other hand, as 
seen from Fig. 4a, the BV5 bulk band changes its position sensitively 
with the photon energy. It reaches the Fermi level near photon ener- 
gies of 21.2 eV (Fig. 4a 1) and 97 eV (Fig. 4a4) where Dirac-like linear 
bands are formed. For other photon energies, the cross point of the 
band appears below the Fermi level. Although we do not know the 
exact K z values corresponding to these photon energies, the obser- 
vation of the Dirac-like bands at some discrete photon energies (K z s) 
is consistent with the theoretical expectations that the conduction 
and valence bands touch only at discrete points in the momentum 
space. As the bulk Cd 3 As 2 has 14 led symmetry identified for the X- 
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Figure 2 | Fermi surface of the Cd 3 As 2 (112) surface and its comparison with band structure calculations, (a). Constant energy contours of Cd 3 As 2 
measured using 21.2 eV light source at a temperature of 45 K. The hexagons depicted in the images represent the (1 12) surface Brillouin zones with the 
first zone plotted as red. They are calculated with lattice constant a = 4.45 A for the (1 12) crystal surface. The measured momentum space covers three 
zone centers. The small dashed hexagon and circle in (a4) and (a5) are guides to the eye. (b). The constant energy contours of Cd 3 As 2 from the band 
structure calculations. The small cell unit calculation takes into account the inverted band structure and modification of the band occupation according to 
the real crystal case. There are two small ring-like Fermi surface sheets around T. With increasing binding energy, the outer sheet becomes warped 
hexagon-shaped at a binding energy of E b = 600 meV while the inner one keeps its ring-like shape. 



ray diffraction patterns, in principle, the Dirac points are along the 
[001] direction 18 due to the C4 symmetry. On the (112) surface, 
however, the Dirac points is located along the [112] direction. This 
is possible because APPES is mostly sensitive to the surface layers 
that has C3 symmetry as observed in the Fermi surface of this surface 
in Fig. 2. It is quite unique for Cd 3 As 2 (112) surface that the bulk 
Dirac cone, as well as the surface Dirac cone, lies close to the Fermi 



level, indicating its electric neutrality that is ideal for investigating the 
intrinsic non-trivial properties. 

In summary, we report angle- resolved photoemission measure- 
ments and band structure calculations on (112) surface of Cd 3 As 2 . 
We find that the overall measured Fermi surface and band structure 
show a good agreement with the band structure calculations. We 
have revealed the existence of two-dimensional Dirac fermions from 
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Figure 3 | Band structures measured along high symmetry cuts for Cd 3 As 2 (112) surface and their comparison with the band structure calculations. 

(a). Location of two momentum cuts along T-M direction (cut A, red line) and T-K direction (cut B, green line), (b). Band structure of Cd 3 As 2 (112) 
surface measured along high symmetry M-r-M direction (cut A in (a)), (c). Band structure measured along M-K-r-K direction (cut B in (a)). To 
highlight the measured bands more clearly, (b) and (c) are second- derivative images with respect to the energy, (d). Calculated band structure along M-K- 
T-M high symmetry lines (blue lines in (a)). Several bulk valence bands (BV) and a surface band (SV) are labeled. 
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Figure 4 | Band structure of Cd 3 As 2 (112) surface measured at different K z and its band structure evolution in the K x -K y plane, (a). Band structure 
along T-M direction (cut A in (c)) measured using different photon energies at a temperature T = 45 K. Different photon energy corresponds to different 
K z . The inner band highlighted with blue dashed line shows weak K z dependence, while the outer band (marked as white dashed lines) shows an obvious 
change with K z . (b). Band structure along several momentum cuts along T-K direction (cuts Bl to B5 in (c)) around T point in the K x -K r plane measured 
using 21.2 eV photon energy. The observed two sets of bands are fitted by hyperbolic function, shown as cyan dotted line for the inner band and 
white dahsed line for the outer band. The velocities used in the fitting functions are 1.4 eV* A and 3.1 eV* A for the outer band and inner band, 
respectively, (c). Location of the momentum cuts in the Brillouin zone. (d). Calculated band structures at different K y in the K x -K r plane. Here, we use the 
same hyperbola functions for the calculated bands as those from fitting the measured bands in (b). 



topological surface state at the (112) crystal surface of Cd 3 As 2 . We 
also present experimental results of 3D Dirac fermions in Cd 3 As 2 by 
observing the linear- dispersive bulk valence band along all the three 
directions in the momentum space. These results are consistent with 
the band structure calculations which indicate that Cd 3 As 2 is a three- 
dimensional Dirac semimetal. Further high resolution ARPES work 
needs to be done to search for the arc-type Fermi surface as expected 
from theory 18 . 
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